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Recently, we outlined a scheme to investigate the effects of compositional order on the magne- 
tocrystalline anisotropy of alloys from a first-principles electronic structure point of view {Phys. 
Rev. Lett. 83, 5369 (1999)} and showed that compositional order enhances the magnitude of mag- 
netocrystalline anisotropy energy (MAE) of Coo.sPto.s alloy by some two orders of magnitude as 
well as affecting the equilibrium magnetization direction. Here we describe our scheme in detail and 
present an in-depth study of the effect by demonstrating its Fermi surface origin. In Coo.25Pto.75 
alloy we find that the perfect LI2 structure has a very small MAE whereas imposition of directional 
order enhances the MAE by two orders of magnitude. We also present the effect of lattice distortion 
(tetragonalization) on the MAE on the same footing and find that in the Coo.5Pto.5 alloy it accounts 
for only about 20% of the observed MAE, thus confirming that compositional order is the major 
player in the enhancement of MAE. We also examine the directional chemical order that can be 
produced by magnetic annealing within the same framework. We extract a Fermi surface mechanism 
for the effect in an explicit study of permalloy. Finally, we propose that the Fermi surface plays a 
major role in the strong coupling between magnetocrystalline anisotropy and compositional order 
in many magnetic alloys. 
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I. INTRODUCTION 

In recent years, there has been intensive theoretical as 
well as experimental research on the magnetocrystalline 
anisotropy of ferromagnetic materials containing tran- 
sition metals, in particular, in the form of multilayers 
and thin films, because of the technological impJieations 
for high-density magneto-optical storage mediaJja Areal 
densities of information storage systems are increasing 
continuously and are expected to reach 40 Gbits per 
square inch by the year 2004 which requires the grain size 
to be less than 10 nm. For such applications the films and 
multilayers need to exhibit a very strong perpendicular 
magnetic anisotropj^'B (PMA) to avoid destabilization 
of the magnetization of recording-bits by thermal fluc- 
tuations and demagnetizing fields. Q Whereas in ultcathin 
films and_miiltilayers the PMA is due to surfaceQ'EI and 
interfacellJllj effects respectively, in thicker films of tran- 
sition metal alloys it is the strong intrinsic bulk mag- 
netocrystalline anisotropy which leads to PMA. Thick 
films of transition metal alloys are particularly interest- 
ing for magneto-optic recording because in addition to 
the required magneto-optic properties these are chemi- 
cally stable and easy to manufacture. To design magnetic 
materials for future magneto-optic recording applications 
a detailed understanding of the mechanism of magne- 
tocrystalline anisotropy is needed. A significant effort 
has been directed towards an understanding of the mi- 



croscopic mechanism of magnetocrystalline anisotropy of 
these systems from a first-principles electronic structure 
point of view but since magnetocrystalline anisotropy 
arises from spin^rbit coupling, which is essentially a rela- 
tivistic effectll3'E3 this means that a fully relativistic elec- 
tronic structure framework is desirable. 

Several experimental observations indicate a correla- 
tion between the compositional order and the magne- 
tocrystalline-anisotropy of ferromagnetic alloys, both m 
the bulkOtll as well as in films.c3 a Very recently,L3 
we developed a 'first-principles' theory of the interrela- 
tionship between the magnetocrystalline anisotropy and 
compositional order, both short- and long-ranged. In 
this theory the electronic structure is treated within the 
spin-polarized fully relativistic Kopringa-Kohn-Rostoker 
coherent-potential approximationEB (SPR-KKR-CPA) 
and the compositional order is modeled using the frame- 
work of static concentralion wavesEll and is an exten- 
sion of our earlier worksE3cj on the magnetocrystalline 
anisotropy of disordered alloys. Within this scheme we 
showed for the first time that in fcc-Coo.sPto.s alloy com- 
positional order indeed has a profound influence on the 
magnetocrystalline anisotropy, especially in the way it 
breaks cubic symmetry. In this paper, we describe our 
scheme in some detail and provide an in-depth study 
of this effect in another CoPt alloy and demonstrate an 
electronic origin of the enhancement of the magnetocrys- 
talline anisotropy energy (MAE). Also, it is known exper- 
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imentally that, upon ordering, the equiatomic CoPt alloy 
undergoes a modest tetragonal lattice distortion (c/a — 
0.98) which can also enhance the MAE. We have cal- 
culated the MAE of disordered fcc-Coo.sPto.s alloy for 
different c/a ratios and find that a 2% tetragonalization 
would contribute only about 20% of the observed MAE. 
This is a clear indication that in CoPt system, the en- 
hancement of MAE is primarily due to compositional or- 
der. This inference gets further credence from the fact 
that fcc-Coo.25Pto.75 alloy, which retains its cubic lattice 
structure upon ordering, also shows an enhancement of 
MAE when compositional order is introduced. 

We have studied the effects of compositional mod- 
ulation on the MAE of fcc-CocPti_c for c=0.25 and 
0.5 alloys. Thick films of these alloys are potential 
magneto-optical recording materials henause they exhibit 
large perpendicular anisotropyj£j"E3o large magneto- 
optic Kerr effect signals compared to the Co/Pt multilay- 
ers and the curceatly used TbFeCo films for a large range 
of wavelength^'El (400-820 nm), have suitable Curie 
temperaturealj (~700 K), high oxidation and corrosion 
resistance as well as being chemically stable and easy to 
manufacture. We investigate the observed ordered struc- 
tures as well as some hitherto unfabricated structures 
of these alloys with same stoichiometry. This is per- 
tinent now that it is possible to tailor compositionally 
modulated films to ahtain better magneto-optic record- 
ing characteristics .CJO We find that, compositional or- 
dering enhances the size of MAE by some two orders of 
magnitude for all the compositions in addition to alter- 
ing the magnetic easy axis in some cases. The easy axis 
in all the layered-ordered structures lies perpendicular 
to the layer-stacking. Of particular interest is the case 
of Coo.25Pto.75. MAE for the perfect LI2 ordered alloy 
is very small as in the homogeneously disordered alloy 
with the easy axis along the [111] direction. However, 
breaking this symmetry by some kind of a directional 
chemical order and creating internal interfaces enhances 
the MAE by two orders of magnitude as well as making 
the easy axis perpendicular to the interfaces, in excel- 
lent agreement with the experimental observations. l3 By 
analyzing the electronic structure of these alloys we find 
that the electrons near the Fermi surface, in particular, 
those electrons around the X-points in the Brillouin zone, 
are largely responsible for the enhancement of MAE by 
compositional ordering. 

Using the same theoretical framework, we have also 
studied the phenomenon of magnetic annealing in 
Nio.75Feo.25 permalloy, which develops uniaxial naagnetic 
anisotropy when annealed in a magnetic field.L3Ed Be- 
cause of its high permeability and low coercivity, the 
permalloy is a good soft magnet and can be used for 
low switching fields in sensors. Also, a large diffeijence in 
the conductivity of majority and minority spinsuj makes 
it a good candidate for spin- valves, spin transistors and 
magnetic tunnel junctions. On magnetic annealing, the 
permalloy develops directional chemical orderCJCJ which 
is responsible for the uniaxial anisotropy. In the previous 



worl£3 we showed for the first time from first-principles 
theoretical calculations that magnetic annealing can in- 
deed produce directional chemical order. Here we isolate 
the electronic origin of the effect. 

The outline of the paper is as follows. In Sec. II, we 
present the formulation. In Sec. Ill we provide some 
numerical and technical details. Afterwards, we present 
our results in Sec. IV and in Sec. V we investigate the 
electronic origin of the MAE enhancement as well as of 
the magnetic annealing. Finally, in Sec. VI we draw 
some conclusions. 



II. FORMULATION 

Our theory is based on the |t:alativistic spin-polarized 
local density functional theopyll'ca and its solution by 
the SPR-KKR-CPA method.Ea'E2l A detailed description 
of the method for the homogeneously disordered ferro- 
magnetic alloys is provided in Ref. pq, and will not be 
repeated. Here, we describe in detail our new frame- 
work to investigate the effects of compositional order, 
both short- and long-ranged, on the magnetocrystalline 
anisotropy (A summary of the approach has already been 
described in our recent letter Ref. |3^). In Sec. II A we 
give a brief outline of the theory of compositional order 
within the SPR-KKR-CPA formahsm and in Sec. II B 
we describe the scheme for studying the effects of com- 
positional order on the MAE together with the theory 
of magnetic annealing. Also, in this section we describe 
a scheme to investigate the electronic origin of the en- 
hancement of MAE upon ordering. 

A. Compositional order 

We consider a binary alloy AcBi^c where the atoms 
are arranged on a fairly regular array of lattice sites. At 
high temperatures the alloy is homogeneously disordered 
and each site is occupied by an A- or i3-type atom with 
probabilities c and (1 — c) respectively. Below some tran- 
sition temperature, Tc, the system will either order or 
phase separate. A compositionally modulated alloy can 
be described by a set of site-occupation variables 
with S^i — 1(0) when the i-th site in the lattice is occu- 
pied by an A(B)-type atom. The thermodynamic aver- 
age, (^i), of the site-occupation variable is the concen- 
tration Ci of an A-type atom at that site. At high tem- 
peratures where the alloy is homogeneously disordered, 
Ci = c for all sites. When inhomogeneity sets in below 
Tc, the temperature-dependent inhomogeneous concen- 
tration fluctuations {6ci} = {ci — c} can bcr-written as a 
superposition of static concentration waves J23 i.e., 
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where, Cq are the ampUtudes of the concentration waves 
with wave-vectors q, and Ri are the lattice positions. 
Usually only a few concentration waves are needed to 
describe a particular ordered structure. For example, 
the CuAu-hke Llg layered-ordered structure (Fig. ^ is 
set up by a single concentration wave with Cq = i and 
q = (001), the [lll]-layered CuPt-like Lli ordered struc- 
ture is set up by a concentration wave with Cq = ^ and 
q = (5^5); and the CusAu-hke LI2 ordered structure is 
set up by three concentration waves of identical ampli- 
tude Cq = J and wave-vectors qi ~ (100), q2 ~ (010), 
and qs = (001) (q is in units of a being the lattice 
parameter). 

The grand-potential for the interacting electrons in an 
inhomogeneous alloy with composition {ci} and magne- 
tized along pLhe direction e at a finite temperature T is 
given by,ErE3 

00 

n{{c,}; e) = i^Z- I de f{e, v) iV({c,}, e; e) 



where, the Onsager cavity correction Ac is given by| 



—00 

^Dc{{ci\;e), 



(2.1) 



where, v is the chemical potential, Z is the total valence 
charge, /(e, v) is the Fermi factor, N{{ci], e; e) is the in- 
tegrated electronic density of states, and il.Dc{{ci}; 
the 'double-counting' correction to the grand-potential.E3 
The derivatives of the grand potential with respect to the 
concentration variables give rise to a hierarchy of direct 
correlation functions. In particular, the second derivative 
evaluated at the equilibrium concentrations. 



c.(2) 
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is the Ornstein-Zecpike direct correlation function for 
our lattice modcETO (so called by the way of the 
close anakjffijf with similar quantities defined for classi- 
cal fluidsE^'EJ) . These are|-£elated to the linear response 
functions, aij{e), throughE^ 



a,y(e) 



/3c(l - c) 



{e)akj{t 



(2.2) 



where /3 — (ksT)^^, fc^ being the Boltzmann constant. 
The linear response functions, aij(e), describe the result- 
ing concentration fluctuations, {Sci}, which are produced 
when a small inhomogeneous set of external chemical po- 
tentials, {Si^i}, is applied at all sites. Via the fluctua- 
tion dissipation theorem these are proportional to atomic 
pair-correlation functions, i.e. aij — f3[{^i^j) — (^i) i d)]) - 
Upon taking the lattice Fourier trarLs£p|tt||i of Eq. ( |2.2| ) 
we obtain a closed form of equations J3 



a{q,T;e) 



/3c(l 



l-/3c(l-c)[5(2)(q;e)-A,]^ 



(2.3) 
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(3c{l - c) Vbz 



— f dq' 5(2)(q'.e)a(q',T;e). 

BZ J 



Here, a{q,T), the lattice Fourier transform of a^j, are 
the Warren-Cowley atomic short-range order (ASRO) pa- 
rameters in the disordered phase. The Onsager cavity 
correction in Eq. (2.3) ensures that th£_sp£ctral weight 
over the Brillouin zone is conserved, E3e30 so that, in 
other words, the diagonal part of the fluctuation dissipa- 
tion theorem is honored, i.e. an = /3c(l — c). 

The spinodal transition temperature Tc below which 
the alloy orders into a structure characterized by 
the concentration wave-vector qmax is determined by 
S'^^^qmax] e), where qmax is the value at which S'^^^q; e) 
is maxim£iL| Ifeglecting the Onsager cavity correction, we 
can write £2I'E3 



c(l-c)5(2)(q„,,;e) 



For the purpose of this study, we need only to consider 
the effects of compositional fluctuations in which charge 
rearrangement effects are neglected. Within the SPR- 
KKR-CPA scheme, a formula for SjiJ{e) is obtained by 
using Ll oyd formulae^ for the integrated density of states 
in Eq. (Q, 



Sfk\^)= -— I def{e,v)Tr {X\e)-X^{e)] 
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where, 

X^(^)(e) 
and 
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X ^ r^"(e) A'"'=(e)T"J(e)j , (2.4) 



0A(i3)(e)-^c-'(e)}"'+r™(e) 



- X'^ie) r^"(e) A"'^(e) T"J'(e) X^{e). 

Here t^(5-)(e) and tc{e) are the t-matrices for electronic 
scattering from sites occupied by A{B) atoms and CPA 
effective potentials respectively, and r'"-' (e) are the path 
operator matrices for the CPA effective medium in real 
space obtained by a lattice Fourier transform of T(k; e), 
where 

r(k;e)= [t-i(e)-5(k)]-\ 

gCk) being the KKR structure constants m atrix .jnl Now 
taking the lattice Fourier transform of Eq. (2.4), we get 
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de fie,v) 



X iL2Ls;LiLi (q; e) AL3L4 (q; e) 

where, 

AL,L,(q;e) = {X^(e)-X^(e)}^^^^ 



(2.5) 



LsLaLsLr 



and, 



3 (q; e) = J dk tl^^l^ (k + q; e) tl^l^ (k; e) 
-^™L3(e)rE^L,(e). (2.6) 

Experimentally, the instability of the disordered phase to 
ordering can be observed in diffuse electron, x-ray, or neu- 
tron scattering experiments. The experimentally mea- 
sured intensities are proportional to the ASRO param- 
eter. In previous works, the ASRO parameter Q;(q, T) 
and ordering temperature Tc have been calculated for 
many alloys, both non-magnetic and ferromagnetic, and 
were pppipared with diffuse x-ray and neutron scattering 
data.E3'E3 However, in those studies, relativistic effects 
were largely ignored. These studies have revealed that 
the electronic structure around the Fermi level is some- 
times the driving force-bpiind unusual compositional or- 
dering in some alloys.EilE^I Schematically, S^^\ci;e) can 
be written in terpis of the convolution of Bloch spectral 
density functionsE3 ^^(k, e; e) as,Eil 



5(2) (q;e) - de de 



f{e,v)- f{e',v) 



-00 —00 



X J dkABik + (i,e;e)AB{k,e';e). (2.7) 

For perfectly ordered alloys, this reduces to a standard 
susceptibility form.E3 In some cases the principal con- 
tribution to 5''^-' (q; e) comes from energies close to the 
Fermi level v. Thus S^^\c[;e) can be large around the 
Fermi energy in two ways.. In, the conventional Fermi 
surface 'nesting' mechanismEilEj overlap takes place over 
extended regions of the reciprocal soace between almost 
parallel flat sheets of Fermi surfaced with spanning vec- 
tor Q as in Cuo.75Pdo.25- The structures based on 
this mechanism will tend to be long period or incom- 
mensurate structures. Alternatively, a large 5'*^^-' (q; e) 
can result from a spanning vector connecting the van 
Hove singularities around high-symmetry points,Ell as in 



Cuo.5Pto.5. The structures based on this mechanism will 
tend to produce high symmetry structures with short pe- 
riodicities since the spanning vector will be that which 
connects the high-symmetry points in the Brillouin zone. 
We find similar Fermi surface effects contribute to the 
enhancement of MAE by compositional order found in 
CocPti_c alloys, and propose that the mechanism is 
widespread in other magnetic transition-metal alloys. 



B. Magnetocrystalline Anisotropy 

The MAE of the inhomogeneous alloy can be charac- 
terized by the change in the electronic grand-potential 
arising from the change in the magnetization direction. 
Thus, 

K{{cA) = mcrhei) - n{{cA;e2), 

where, ei and 62 are two magnetization directions. We 
assume that the double-counting correction ^Dci{ci}, e) 
is generally unaffected by the change in the magnetiza- 
tion direction, and therefore, only the first two terms of 
Eq. (O) contribute to the MAE, 



Ki{c^}) = [lyi - V2)Z - J de fie, vi) 7V({c,}, £; ei) 

— 00 

00 

+ / de f{e,v2) N{{cA,e-e2), 



where vi and 1^2 are the chemical potentials of the sys- 
tem when the magnetization is along ei and 62 directions 
respectively. The change in the chemical potential origi- 
nates from a redistribution of the occupied energy bands 
in the Brillouin zone in the event of a change of magne- 
tization direction. A Taylor expansion of /(e, V2) about 
vi and some algebra leads to. 



K{{cA)^ / d£/(£,i.i)[7V({c,},£;ei)-A^({cJ,£;e2)] 



+ 0(i/i - V2f, 

Note that the effect of the small change in the chemi- 
cal potential on K{{ci}) is of second order in [vi — 1^2), 
and cajji be shown to be very small compared to the first 
term.E3 We now expand K{{ci}) around Kcpa{c), the 
MAE of the homogeneously disordered alloy AcBi-c, 



K{{cA)=Kcpa{c)+J2 



dK{{c,}) 



dc 



6ci 



IE 



d'mcA) 



2 ^ dcjdck 



5cj5ck + 0{5cf . (2., 
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It is clear from Eq. (2.8) that we are considering only the 
effects of two-site correlations on the electronic grand- 
potential and MAE. However, in principle it is possible 
to include higher order correlations, but computationally 
it will be prohibitive. 

Within the SPR-KKR-CPA scheme, a formula, for 
Kcpa{c) is obtained by using the Lloyd formulaic for 
the integrated density of states, 

oo 

— OO 

xln||/+{t-i(e2)-<-\ei)}r(k;ei)|| 
+ c (lnP^(ei)||-ln||Z?^(e2)||) 

+ (l-c) [ln\\D^ie,)\\-\n\\D^{e2)\\)'_ 
+ 0{yi-V2f, (2.9) 



where 

i?-4(S)(e) = 



rO°(e){t-i 



MB) 



(e)-ie-'(e)} 



Note that Eq. (2_^) is the finite temperature version of 
the expression of MAE of disordered alloys given in Ref. 
|38| . Now, the variation of the MAE with respect to the 
change in the concentration variables is. 



dK{{c.}) 



da 



JjYl f 

- de f{e,vi) 

— OO 

lnp^(ei)||-lnp^(ei) 
An\\D^{B2)\\+\n\\D^{e^) 



which is indep end ent of the site-index and so the second 
term in Eq. (2^) vanishes if the number of A and B 
atoms in the alloy is to be conserved ( dcj = ). 
Also, we have. 



dcjdck 



(2) 

where, S jjj {sk^-t^ are the Ornstein-Zernike direct corre- 
lation functionsLj when the magnetization i s alo ng ei(2)- 
Now taking the Fourier transform of Eq. (2.8), we get 
the MAE of the compositionally modulated alloy with 
wave- vector q. 



K{ci) - Kcpa{c) 



1 



i^(2)(q). 



where 



(2.10) 



(2.11) 



Eq. ( ^.lOD thus shows a direct relationship between the 
type of compositional modulation and the MAE. Also, by 



calculating S^'^^ (q; e) for different q- vectors, while keep- 
ing the magnetic field and magnetization direction fixed, 
one can study the effect of an applied magnetic field on 
the compositional modulation of a solid solution, and 
thus can describe the phenomenon of magnetic annealing. 
In this case, the q- vector for which 5*^^-' (q; e) is maximal 
will represent the compositional modulation induced in 
the alloy when it is annealed in the magnetic fjeM. . 

In our previous studies on disordered alloya23c3 we 
have found that the MAE originates from the change in 
the electronic structure around the Fermi energy caused 
by altering the magnetization direction. The changed 
electronic structure around the Fermi surface caused by 
the compositional order is thus the reason for the en- 
hancement of MAE in the ordered phase. To understand 
this effect we undertake the foll owin g analysis. Using the 
same argument leading to Eq. (2/7) we arrive at, 



OO OO 



i4:(2)(q) _ J de J de 



e — e' 



— OQ — OO 



X J dk [AB(k + q,e;ei) AB(k,e';ei) 
-^^(k + q, e; 62) ^^(k, e'; 62)] 

which can be rewritten as, 

00 00 



e - e' 



—00 —00 



where, 



X J dk [S(k-|-q,e) A(k,e') 

+ E(k,£') A(k + q,e)] (2.12) 

S(k,£) = AB(k,£;ei)+AB(k,e;e2) (2.13) 



and 



A(k,£) :=AB(k,e;ei)-^B(k,e;e2). 



(2.14) 



The principal contributions to K^^^ (q) will come from en- 
ergies close to the Fermi energy v. Thus, i^(2)(q) will be 
large for q- vectors corresponding to the spanning vectors 
q = Q connecting the peaks in I](k, e) and A(k, e) for 
e close to the Fermi energy. Therefore, a close examina- 
tion of S(k, v) and A(k, will give an idea as to which 
composition modulations will produce a large MAE. 



III. COMPUTATIONAL DETAILS 

In this section we discuss the technical details in the 
calculation of MAE of compositionally modulated alloys. 
In a homogeneously disordered alloy, the MAE is of the 
order of ^eV and the MAE of the ordered systems can 
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be as high as 0.5meV. But the MAE is still several or- 
ders of magnitude smaller than the band energies as well 
as the total energy. For this reason the calculation of 
MAE needs great care. In an earlier publicatioiiE3 we 
presented a robust scheme to calculate the MAE, where 
we proposed that one should calculate the difference in 
the total energies for the two magnetization directions 
directly rather than subtracting the two total energies 
calculated separately. We included a detailed discussion 
on the computational and technical aspects of solving 
the SPR-KKR-CPA equations as well as the calculation 
of MAE of the disordered alloys and do not repeat that 
discussion here. Instead, we will focus only on the calcu- 
lation of S'(2)(q;e) and K'^^\ci). 



T he int egration over the energies involved in Eqs. (2.5) 
and (2.11) is carried out using a complex contour because 
in the complex plane, the integrand becomes a smooth 
function of both energy as well as k and consequently 
we need fewer energy points-as well as fewer k-points to 
obtain an accurate integral.t3 In our calculation, we have 
used a rectangular box contour as described in detail in 
Ref. ||. 

The most demanding part of the whole calculation is 
t he evaluation of the convolution /(q; e) given by Eq. 
(2.6). Owing to the form of the integrand the integra- 
tion has to be done using the full Brillouin zone. We 
performed the Brillojuin zone integration by using the 
adaptive grid methodEj which has been found to be very 
efficient and accurate. In this method one can preset the 
level of accuracy of the integration by supplying a toler- 
ance parameter e. The number of k-points used in the 
integration thus depends directly on e. In the present 
work, the integration in Eq. ( |2.6| ) is carried out with 
e — 10~^ which means that values of S'*^^-'(q; e) which are 
of the order of 0.1 eV are accurate up to 0.1 ^eV. To 
achieve such level of accuracy, we had to sample a large 
number of k-points in the Brillouin zone. Typically, in 
our calculations, we needed around 35,000 k-points in the 
full Brillouin zone for energies with an imaginary part of 
0.5 Ry. When the energy was 0.0001 Ry above the real 
axis (and that is the closest point to the real axis, we 
need) we required as many as 3 million k-points in the 
full Brillouin zone for the same level of accuracy. Fur- 
thermore, we have calculated S''^^-'(q; ei) and 5''^-'(q;e2) 
simultaneously ensuring that they are calculated on the 
same grid of k-points, and thus cancelling out the sys- 
tematic errors if there are any. Therefore, we claim that 
the values of K{q) are accurate to within 0.1 /ieV. 



IV. RESULTS AND DISCUSSION 

A. CocPti_c alloy 

We have studied CocPti_c alloys for c=0.5 and 0.25. 
It is well known that disordered fcc-Coo.sPto.s under- 
goes a phase transformation into a CuAu-type LIq or- 



dered tetragonal structurelliS with a c/a ratio of 0.98 
below 1100 K, and Co(i:^Pto.75 orders into a CuaAu- 
type LI2 cubic structural below 960 K. Thus, in case of 
the equiatomic composition, the tetragonalization of the 
lattice which also lowers the symmetry could also be re- 
sponsible for MAE enhancement, whereas in Coo.25Pto.75 
there is no such additional effect, because the lattice re- 
mains cubic even in the ordered phase. With this in mind 
we calculated the MAE of disordered volume-conserving 
face-centered-tetragonal Coo.5Pto.5 alloy as a function of 
the c/a ratio. The results are presented in Fig. ||. We 
point out that in these calculations we have used atomic- 
sphere approximation for the single-site potentials and 
that the potentials and the Fermi energy are those of 
the disordered fcc-Coo.5Pto.5 alloy. To estimate the mag- 
nitude of the effect of tetragonalization we have 'frozen' 
these potentials as the c/a ratio has been altered. We ob- 
serve that the MAE is a monotonically decreasing func- 
tion of the c/a ratio, and that it is positive for c/a < 1.0 
and negative for c/a > 1.0. This is consistent with the 
experimental observations that the magnetostriction con- 
stant, Aooi, is positiveEil'E3 foir4he disordered Coo.5Pto.5 
alloy. Because, Freeman et aBH have shown that Aooi is 
proportional to the rate of change of MAE with respect to 
the c/a ratio, and that Aqoi has the opposite sign to that 
of the later (Note that, there is a sign difference in our 
definition of MAE and that of Freeman et at). There- 
fore, our results are in qualitative agreement with the 
experimental observations. The MAE at the experimen- 
tal value of c/a (0.98) is about 25 ^eV which is less than 
20% of the experimentally observed MAE. Most impor- 
tantly, the sign of MAE for this value of c/a is positive 
which means that the magnetic easy axis is not along 
the [001] direction (c-axis) in direct contradiction to the 
experimental observations. Consequently, we conclude 
that in the ordered alloys of Co and Pt lattice distortion 
is not the major factor in determining the MAE rather it 
is the compositional order which is primarily responsible 
for the large MAE. 

Now we present the effect of compositional order on 
the magnetocrystalline anisotropy. In our studies, we ex- 
plore the observed equilibrium ordered structures as well 
as some hypothetical ordered structures of these alloys 
keeping the stoichiometry constant. The results are sum- 
marized in Tables | and ||. First we discuss Coo.sPto.s. 
We note that, S''^^)(q) is maximum for the LIq structure, 
implying that the alloy would order into this structure 
at 1360 K as it is cooled from high temperature in good 
agreement with experiment (experimental value of the 
ordering temperature is 1100 K). From Table || we note 
that for q = (001) and q — (100) which correspond to 
CuAu-like LIq ordered structure, with Co and Pt layers 
stacked along the [001] and [100] directions respectively, 
the direction of spontaneous magnetization is along the 
[001] and [100] directions respectively in excellent agree- 
ment with experiment. Also, the MAJiL(S8.6 /ieV) is com- 
parable to the experimental valueliJiij (~ 130/ieV) as 
well as to the calculated value of the Llo-ordered tetrag- 
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onal CoPt alloy.Ea In the CuPt-type Lli layered structure 
set up by q = (555), the equihbrium duection of mag- 
netization is along the [111] direction of the crystal, i.e. 
again perpendicular to the layer-stacking. The structure 
set up by concentration waves with q — (10^) and (01^) 
is also a [001]-oriented layered structure, but the layers 
are not alternately pure Co and Pt planes, rather they are 
layers of CoPt in a particular order (Fig. |l|). Even in this 
case, we note that the magnetization is perpendicular to 
the layered structure. Similarly, for q — (5OI) and (5IO) 
where the planes are stacked along the [100] direction the 
magnetization is also along the [100] direction. 

In case of Coo.25Pto.75, 'S'^^^(q) has a maxima for 
q=(100), (010), and (001) implying that a LI2 ordered 
structure is favored with a transition temperature of 935 
K in excellent agreement with the experimental value of 
960 K (Ref. ^5|). Note that, a combination of three wave- 
vectors qi = (100), <l2 ~ (010), and qs — (001) generates 
the isotropic LI 2 ordering, while a single wave- vector, 
for example, qi = (001) generates a layered structure 
with directional compositional ordering along the [001] 
direction. This is a superstructure consisting of alternat- 
ing monolayers of pure Pt and Coo.5Pto.5, as depicted 
in Fig. ^ In this structure, therefore, there are no 
out-of-plane Co-Co bonds, only in-plane Co-Co bonds 
which can produce only in-plane Co-Co nearest neighbor 
pairs. We find that for this structure, the easy axis is 
along the [001] direction, i.e. perpendicular to the layers, 
and the MAE is also quite large (~ 84/ieV). In a recent 
experimenttJ it was found that, the [001]-textured thick 
films of Coo.25Pto.75 alloy deposited at 670 K do have 
this type of structure, i.e. there are stacks of Pt and 
Coo.5Pto.5 monolayers perpendicular to the [001] direc- 
tion, and these films exhibit PMA. It should be empha- 
sized, however, for the perfect LI2 structure, where all 
the three wave-vectors, namely, qi = (100), q2 = (010), 
and q3 = (001) contribute, the easy axis is along the 
[111] direction and the MAE is very small, comparable 
to that of the disordered alloy. Another [lll]-stacked 
layered structure is generated by q — (55^) in which 
there are alternating planes of pure Pt and Coo.sPto.s. 
In this case also the easy axis is perpendicular to the 
layer stacking (i.e. along [111]) and the MAE is 23.4 
/ieV. We are not aware of any experimental results on the 
bulk ordered Coo.25Pto.75 system. However, it is reported 
that [lll]-textured thick films grown around 690 K hav- 
ing anisotropic .compositional order exhibit large uniax- 
ial anisotropycJHS whereas films deposited around 800 
K with a Ll2-type isotropic chemical order exhibit no 
anisotropy. These observations are clearly in good agree- 
ment with our results. 

The above observations suggest that a large MAE and 
an easy axis perpendicular to the layer-stacking is a result 
of existence of in-plane Co-Co bonds and out-of-plane 
Co-Pt bonds. It has been observed experimentally that 
in these films there are indeed more Co-Co bonds in- 
plane and Co-Pt bonds out-of-plane which produce in- 
ternal interfaces analogous to Co/Pt multilayers. E3 Be- 



cause of the out-of-plane orientation of Co-Pt pairs a hy- 
bridization between the Co atoms with large magnetic 
moment and Pt atoms with strong spin-orbit coupling in- 
duces anisotropics in the 3d and 5d orbital moments, the 
out-of-plane components of the orbital moments become 
larger than the in-plane components, as has been ob- 
served in the x-ray magnetiajpircular dichroism measure- 
ments in the CoPts films.E3o The induced anisotropy in 
the orbital moments directs the spin moment into a per- 
pendicular alignment through spin-orbit coupling, thus 
overcoming the in-plane-shape anisotropy due to the spin- 
spin dipole interactionEa 



B. Magnetic annealing of Nio.75Feo.25 alloy 

When Nio.75Feo.25 permalloy is annealed in a mag- 
netic field, an uniaxial magnetic anisotropy is induced 
depending on the direction of the apiplicd field with re- 
spect to the crystallographic axes.EErEil This phenomenon, 
known as magnetic annealing, is interpreted as the. cr.e-| 
ation of directional chemical order in the material Jlj'c£l 
A recent stujdy based on magneto-optic Kerr effect 
measurementalJ also reveals that magnetic annealing can 
induce uniaxial anisotropy. Previous studies based on 
non-relativistic electronie-, structure calculationsO and 
Monte Carlo simulationgl3 have shown that in this sys- 
tem compositional and magnetic ordering have a large 
influence on each other. However, to describe magnetic 
annealing a fully relativistic treatment is needed. We 
produce the first quantitative description of magnetic an- 
nealing from ab initio electronic structure calculations in 
Nio.75Feo.25- Our results are summarized in Table IH. 
We calculate S'(^^(q, e) for the permalloy in an appHed 
magnetic field of strength 600 Oe (same as used in the 
experimented) oriented along e = [001], [111], and [100] 
directions of the crystal. We find that when the mag- 
netic field is along [001] (column 2) 5^^^ (q) is maximum 
for q=(001) confirming that ordering is favored along 
the direction of applied field. In this case, the compo- 
sitional structure is a layered structure along the [001] 
direction comprising of alternate pure Ni layers and dis- 
ordered Nio.5Feo.5 layers. Noting that the measured in- 
tensity in a scattering experiment is proportional to the 
ASRO parameter a{q), we present the ASRO parame- 
ters calculated at 722 K (1 K above the ordering tem- 
perature) for different q- vectors in Table IV. We note 
that when the applied magnetic field is along the [001] 
direction, a(q) for q=(001) is 40-50% larger than that 
at other q-vectors. Therefore, when the alloy is annealed 
in the magnetic field, the superlattice spot in the mea- 
sured intensity at q=(001) will be 40-50% more intense 
than that at q=(100) at a temperature 1 K above the 
transition temperature. When the magnetic field is ap- 
plied along the [100] direction, a similar structure per- 
pendicular to [100] direction is favored. However, when 
the applied field is along [111] direction (column 6) all 
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the three orderings, namely, (100), (010), and (001) are 
favored equally. Thus, in this case, we will get a CuaAu- 
type LI2 ordering. The calculated transition tempera- 
ture 221 K is in good agreement with the experimental 
valued of 820 K. 

We have also calculated the MAE of Nio.75Feo.25 
permalloy for different ordered structures. We found that 
the MAE of the disordered phase is very small (less than 
0.1/xeV) and in case of directional chemical ordering it 
is of the order of /zeV, an increase by an order of mag- 
nitude. However, in case of perfect ordering, which is 
LI2, the MAE becomes very small as in the disordered 
phase. This is-consistent with [experimental observations 
of ChikazumiEj and Ferguson. 



V. ELECTRONIC ORIGIN OF MAGNETIC 
ANISOTROPY 

In elemental solids and the disordered alloys the 
MAE originates from a redistribution of electronic states 
around the Fermi lei2B|l caused by the change in the mag- 
netization direction.Ea The electronic structure of the dis- 
ordered phase around the Fermi level is also partly re- 
sponsi|bie--£or the tendency to compositional order in some 
alloysEiI'LJ The enhancement of MAE in the composi- 
tionally modulated alloys is also related to the electronic 
structure around the Fermi level in the disordered phase. 
Consequently, to understand the electronic origin of large 
MAE in the ordered phase, in the following we analyze 
the Bloch spectral function of the disordered Coo.5Pto.5 
alloy along the lines as depicted in Eq. (2.12). 

First we consider the ordering tendency. As presented 
in Table ||, our calculation predicts Llo-type order in 
Coo.5Pto.5 which is in excellent agreement with exper- 
iment. The quantity ^(^^(q; e) given by Eq. ( |2.5| ) can be 
rewritten as. 



oc 

5(2)(q.e) = -:^ J de f{e,i,)Fiq,s;e) 

— 00 

where, 

F{<i,e;e)^ [{X^{e) - {e)} ^^^^ 

L 1 Zi2 i'3 i'4 

X /L2L3;L4ii(q;e) AL3L,(q;e) . (5.1) 

S'^^^ (q) is an integiiated quantity and can also be_sKritten 
in terms of a surrO over Matsubara frequencies Jl3 w„ — 
{2n + l)nkBT, 

^(2) (q; e) = 2kBT ^ Re [F(q, + icu^; e)] . 

n 

In Fig. H we show a plot of F(q, e; [001]) calculated for 
energies along the imaginary axis perpendicular to the 



Fermi level for q-vectors (000), (001), (iii) and (10^). 
The area under these curves is indicative of the strength 
of S*'^' (q; [001]) for these q-vectors and is obviously great- 
est for q=(100). 

In special cases the electronic structure near the Fermi 
surface can explain some unusual ordering tenden cies .. 
One famous example is a nesting of the Fermi surfaccHS 
which takes place over extended regions of the reciprocal 
space between almost parallel sheets of Fermi surfaces 
as in Cuo.75Pdo.25 and gives rise to its incommensurate 
ordering tendency. However, presence of van Hove sin- 
gularities can give rise to a somewhat different kind of 
mechanism of ordering in which the spanning vector cop-, 
pies only the regions around two high-symmetry pointsEil 
as in the case of CuPt. In Fig. ||we show a density-plot 
of the Bloch spectral density function of the disordered 
Coo.5Pto.5 alloy in the (001) plane {k^ = 0). In this figure 
the F-points are at the corners and the X-points are at 
the center of the edges. White areas indicate relatively 
larger density of states. The cross-section of the Fermi 
surface sheets are seen quite clearly. In this figure we do 
not observe any van Hove singularities at the symmetry 
points, the ordering tendency coming instead from states 
over a wide energy range. However, we do notice that, 
there are large densities of states around the X-points as 
well as around the (110) points. 

Now we discuss the enhanceme nt of MAE. Again, the 
quantity K^^'>{q) given by Eq. ( 2.11 ) is an integrated 



quantity which can be written as a sumxif contributions 
evaluated at the Matsubara frequencies ,E3 

if (2) (q) ^2kBTj2 Re [AF{(i, v + (5.2) 



where 



Ai^(q, e) = F(q, e- ei) - i^(q, e; 62). 



(5.3) 



In Fig. H we show a plot of AF(q, e) with ei = [001] 
and 62 — [111] calculated for energies along the imagi- 
nary axis perpendicular to the Fermi level for q-vectors 



(000), (001), ( 



111^ 

2 2 2> 



and (lOi). The area under these 



curves is indicative of the magnitude of if (q) for these 
q-vectors. The striking feature is that the principal con- 
tributions are near the real axis indicating that the Fermi 
surface plays the dominant role in the enhancement of 
the MAE as well as the direction of easy magnetiza- 
tion. We propose that such Fermi surface contributions 
are widespread factors in the enhancement of MAE by 
compositional order in many magnetic alloys. 

The origin of the enhancement of MAE and the role 
of the Fermi surface can be easily understood by ana- 
lyzing the sum and difference of Bloch spectral density 
function s for the two magnetization directions as outlined 
in Eq. ( 2.12 ). Obviously, we are looking for the span- 
ning vector Q which connects the high density regions 
in I](q, t^), the sum of the Bloch spectral density func- 
tions for two magnetisation directions, and A{q,v), the 
difference of the Bloch spectral density functions for two 
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magn etisation directions, as defined in Eqs. ( 2.13| ) and 
( ^.14 ). In Fig. ^we show a density-plot of E(q, ^) and 
A(q, i/) of the disordered Coo.sPto.s alloy at the Fermi 
energy in the (010) plane {ky = 0). In this figure the 
F-points are at the corners and the X-points are at the 
center of the edges. In the density-plot of S(q, i^) the 
white areas indicate relatively larger values and the most 
dark areas indicate zero values, whereas in the density- 
plot of A(q, v) the most white (most dark) areas indicate 
large positive (negative) values and the grey areas indi- 
cate intermediate values. As in case of the Bloch spectral 
density function, we do not observe large values at the 
symmetry points. The analysis is more complicated for 
other q-vectors but we can conclude that the enhance- 
ment of MAE for concentration wave-vectors near the 
Brillouin zone boundary is indeed due to large values of 
S(q, i^) and A(q, i^) at the Brillouin zone boundary. 

Now we discuss the electronic origin of magnetic an- 
nealing effect in Nio.75Feo.25 alloy. We note from Ta- 
ble that the difference between S'^^) (q; [001]) and 
^(^^(q; [100]) for different q-vectors is quite small, but 
as depicted in Table IV it is large enough to be observed 



by diffuse x-ray, electron and neutron scattering experi- 
ments. Because of the small difference in these quantities, 
in Fig. we plot the real part of AF(q, e) which is related 
to the difference between ^(^^(q; [001]) and S'^^Hq; [100]) 
for q=(001) for energies along the imaginary axis per- 
pendicular to the Fermi level (left-hand scale). We also 
plot i^(q,e;[001] which is related to S'^^^Q; [001]) for 
q=(001) for the same energies (right-hand scale). We 
note that, these two quantities peak near the Fermi en- 
ergy, i.e. when the imaginary part of the energy is very 
small. This is an indication that this process is also gov- 
erned by the electrons near the Fermi surface. 



[Ill] direction which is the same as in its disordered coun- 
terpart. The MAE is also not very large. However, in 
case of directional ordering along either of the [100], [010] 
or [001] direction the easy axis is along the ordering with 
greatly enhanced MAE. Finally, by analyzing the elec- 
tronic structure of the disordered alloy near the Fermi 
energy we have found that the Fermi surface plays the 
dominant role in the enhancement of MAE. Within the 
same theoretical framework we have also been able to 
explain the appearance of directional chemical order in 
Nio.75Feo.25 when it is annealed in an applied magnetic 
field and linked that also to the alloy's Fermi surface. As 
a last remark, we look forward to future work in which 
the effects of compositional structure are fully incorpo- 
rated into micromagnetic modeling of transition metallic 
materials via ab initio electronic structure calculations. 
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VI. CONCLUSIONS 

We have presented the details of our ab initio theory of 
the connection of magnetocrystalline anisotropy of ferro- 
magnetic alloys with the compositional order within the 
SPR-KKR-CPA scheme. This theory has been applied 
to fcc-CocPti_c alloys. We found that when cooled from 
a high temperature, fcc-Coo.5Pto.5 tends to order into 
ilo layered-ordered structure around 1360 K and fcc- 
Coo.25Pto.75 tends to order into LI2 structure around 935 
K, in good agreement with experimental observations. 
Also, we found that in Llo ordered CoPt the sponta- 
neous magnetization is along the [001] direction which 
assumes the c-axis of the tetragonal structure, in excel- 
lent agreement with experiment. The magnitude of the 
MAE is also close to the experimental value. In some 
of the other hypothetical layered-ordered structures with 
same stoichiometry also the magnetic easy axis lies per- 
pendicular to the layer stacking, with the [lll]-stacked 
CuPt-like Lli structure having the largest MAE. In the 
LI2 ordered CoPta the magnetic easy axis is along the 
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TABLE I. Direct correlation function ^(^'(q; [001]) for dif- 
ferent q-vectors for CocPti_c alloys (the respective ordered 
structures are shown in Fig. hh. 



c 


q 


Structure 


^(^'(q; [001]) 
(eV) 


(K) 




(000) 


Clustering 


-1.51 






(001) 


Llo 


0.47 


1360 


0.5 


(100) 

(---) 

(101) 


Llo 
LU 


0.47 
0.29 
0.19 


1360 




(100) 


LI2 


0.43 


935 




(010) 


LI2 


0.43 


935 


0.25 


(001) 

(111) 

^2 2 2' 

(llo) 
(iOl) 


LI2 
DO22 


0.43 
0.22 
0.19 


935 




DO22 


0.19 





11 



TABLE II. Magnetocrystalline anisotropy energy K{q) for several compositionally modulated CocPti_c alloys characterized 
by different q- vectors (the respective ordered structures are shown in Fig. |l|). Here A'(q) are calculated with respect to the 
reference system which has the magnetization along the [001] direction (i.e. ei — [001]) of the crystal. Thus, when A'(q) < 



the easy 


axis is along [001] and when -K'(q) > the easy axis is along 


: 62. 






c 


q 


Structure 


B2 — [lllj 

Kid) 


^(q) 

(^eV) 


Easy Axis 


0.5 


(001) 
(100) 

(III) 

(lOi) 

\^ 2' 

(|01) 


Llo 

Llo 

LU 
[001]-Layered 
[100]-Layered 


-58.6 
39.6 
152.0 
-158.7 
85.0 


-105.6 
105.9 
0.0 

-236.5 
236.3 


[001] 
[100] 

[111] 

[001] 
[100] 


0.25 


(100) 
(010) 
(001) 

(111) 

(100), (010), (001) 
(liO) 
(Oli) 


[100]-Layered 
[010]-Layered 
[OOlj-Layered 
[lll]-Layered 
LU 


30.5 
30.5 
-53.0 
23.4 
8.0 
-52.4 
100.2 


83.6 
0.0 

-83.6 
0.0 
0.0 
0.0 

152.5 


[100] 
[001] 

[111] 
[111] 

[001] 
[100] 
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TABLE III. Direct correlation function 
ferent q-vectors for Nio.75Feo.25 alloy (the 
structures are shown in Fig. |l|). 



S'<^'(q;e) for dif- 
respective ordered 



e=[001] 
(meV) (K) 



e=[100] 
(meV) (K) 



e=[lll] 
S(')(q) n 
(meV) (K) 



(100) 
(010) 
(001) 

(111) 

V 2 2 2 

(lOi) 
(iOl) 



330.880 
330.880 
331.158 
-99.291 
99.284 
99.096 



720.9 



331.158 
330.880 
330.880 
-99.291 
99.096 
99.285 



720.9 



330.973 
330.973 
330.973 
-99.475 
99.222 
99.222 



720.5 
720.5 
720.5 



TABLE IV. Atomic short-range order (ASRO) parame- 
ter a(q; e) calculated at temperature 1 K above the order- 
ing temperature 721 K at few high-symmetry q-vectors for 
Nio.75Feo.25 alloy (the respective ordered structures are shown 
in Fig. i). 



q 


e=[001] 


e=[100] 


e=[lll] 


(100) 


1509.7 


2601.4 


1756.2 


(010) 


1509.7 


1509.7 


1756.2 


(001) 


2601.4 


1509.7 


1756.2 


(111) 

V 2 2 2 ^ 


2.3 


2.3 


2.3 


(lOi) 


4.3 


4.3 


4.3 


(iOl) 


4.3 


4.3 


4.3 



AB 



Llo:q=(0 1) Llj: q=(i 



q=(10i) 



A 





/ 








® 




® 
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Lh: qi=(l 
q2=(0 1 
q3=(0 1 



q=(0 1) q=(i i i) 



FIG. 1. Some ordered structures and their representative 
concentration wave-vectors. For the A_B-type stoichiometry 
q=(001) and q=(m) generate respectively the CuAu-type 
Llo layered ordered structure with layers perpendicular to 
the [001] direction and the CuPt-type Lli layered-ordered 
structure with layers perpendicular to the [111] direction 
and q=(10^) generates a layered structure with planes of 
an ordered structure of A and B atoms stacked along the 
[001] direction. For the AB^-type composition, a com- 
bination of qi = (100), q2 = (010), and q3 = (001) generates 
the CusAu-type LI2 ordered structure. For this composi- 
tion, a single wave-vector q=(001) generates a superstruc- 
ture of alternating monolayers of pure B atoms and disor- 
dered v4o.5-Bo.5 perpendicular to the [001] direction. Similarly, 
q=(iii) generates a superstructure of monolayers of pure B 
atoms and disordered j4o.5-Bo.5 perpendicular to the [111] di- 
rection. The full circles denote A atoms, open circles denote 
B atoms, and a full circle inscribed by an open circle denotes 
a CPA effective atom y4o.5-Bo.5- 
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FIG. 2. Magnetocrystalline anisotropy energy (MAE) of 
disordered volume-conserving face-centered tetragonal (fct) 
Coo.sPto.s alloy as a function of the c/a ratio. 



FIG. 4. Density-plot of the Bloch spectral density function 
of disordered Coo.sPto.s alloy at the Fermi energy in the (001) 
plane (fcz = 0). The F-points are at the corners and the 
X-points are at the center of the edges. White areas indicate 
relatively larger density of states. 
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FIG. 3. Real part of F(q,e, [001]), as given by Eq. (|5.1), 
for q=(100) (full line), (iii) (long-dashed line), (lo|) 
(short-dashed line), and (000) (dotted line) at complex en- 
ergies along the imaginary axis perpendicular to the Fermi 
level for disordered Coo.sPto.s alloy. 
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FIG. 5. Real part of AF(q,e), as defined in Eq. ( |5.3[ ), 
with ei = [001] and 62 = [111] for q=(100) (full line), (001) 
(long-dashed line), (|^|) (short-dashed line), and (000) (dot- 
ted line) at complex energies along the imaginary axis per- 
pendicular to the Fermi level for Coo.sPto.s alloy. Note the 
different scales on the y-axes of this figure and Fig. 0. 
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FIG. 6. Density-plot of the sum over two different mag- 
netization directions of the Bloch spectral density functions 
E(q, v) (top) and the difference of the Bloch spe ctral den- 
s ity f unctions A(q, i/) (bottom) as defined in Eqs. (2.13) and 



( 2.14 ), for disordered Coo.sPto.s alloy at the Fermi energy in 
the (010) plane {ky — 0). The F-points are at the corners and 
the X-points are at the center of the edges. White areas indi- 
cate relatively larger density of states. In the bottom figure 
the dark areas indicate negative values rather than zero. 
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FIG. 7. Real part of AF{q, e) as defined in Eq. (U) with 
ei = [001] and eg = [100] (full line , lef t scale) and real part 
of F(q,e; [001]), as defined in Eq. (|1]) for q=(001) (dashed 
line, right scale) at complex energies along the imaginary axis 
perpendicular to the Fermi level for Nio.75Feo.25 alloy. 
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